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ABSTRACT 

We use the data of Wisconsin Ha Mapper (WHAM) to test the hypothesis of whether the 
amplitudes and spectrum of density fluctuations measured by WHAM can be matched to the 
data obtained for interstellar scintillations and scattering. To do this, first of all, we adjusted 
the mean level of signal in the adjacent patches of the data. Then, assuming that the spectrum 
is Kolmogorov, we successfully matched the amplitudes of turbulence obtained from the WHAM 
data and the interstellar density fluctuations reported in the existing literature. As a result, we 
conclude that the existing data is consistent with the Kolmogorov cascade which spans from 10 6 
to 10 17 m. 

Subject headings: methods: data analysis — turbulence — ISM: lines and bands — techniques: 
spectroscopic 

1. Introduction 

The interstellar medium (ISM) is turbulent on scales ranging from AUs to kpc (see Elmegreen & Scalo 
2004 for a review), with an embedded magnetic field that influences almost all of its properties. MHD 
turbulence is generally accepted to be of key importance for fundamental astrophysical processes, e.g. star 
formation, propagation and acceleration of cosmic rays (McKee & Ostriker 2007, Lazarian et al. 2009). 

A substantial insight into the scale density fluctuations has been obtained via studies of radio scattering 
and scintillations. It was shown in Armstrong et al. (1995) that the spectrum of density inhomogeneties is 
consistent with the Kolmogorov turbulence over a wide range of scales from 10 6 m to 10 13 m (see also Spangler 
& Gwinn 1990, Wang et al. 2005, You et al. 2007). The theoretical explanation of this effect is provided 
within Goldrcich-Sridhar (1995) model of turbulence, where Kolmogorov spectrum of magnetohydrodynamic 
fluctuations is expected (see also Cho & Lazarian 2003, Beresnyak & Lazarian 2009). But what is happening 
at larger scales? 

The question was addressed by a number of recent studies. For instance, Haverkorn et al. (2008) used 
the Southern Galactic Plane Survey data and concluded that the injection scale of the turbulence in the 
spiral arms of the Galaxy should be rather small, e.g. order of a parsec, while a much larger, e.g. order 
of 100 pc, outer scale is reported for the interarm regions. This, however contradicts to the smooth power 
laws observed for velocity fluctuations in HI (Green 1993, Lazarian & Pogosyan 2000) and thus stimulates a 
further quest of density fluctuations at large scales. 

In the paper below we address this question by analyzing the Wisconsin Ha Mapper (WHAM) data 
(see Madsen, Haffner & Reynolds 2002). We check the consistency of the fluctuations of WHAM and those 
measured by scintillations and scattering to the Kolmogorov cascade. 
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In what follows, we discuss the WHAM data in Sect. 2, expected angular spectrum of fluctuations in 
Sect. 3, determination of the spectrum's free parameters in Sect. 4 and provide discussions in Sect. 5. In 
Appendix A we describe the algorithm of correcting background level errors. 



2. WHAM data and analysis 

The Wisconsin Ha Mapper has surveyed the distribution and kinematics of ionized gas in the Galaxy 
above declination -30°. The WHAM Northern Sky Survey (WHAM-NSS) has an angular resolution of 1° 
and provides the first absolutely calibrated, kinematically resolved map of the Ha emission from the warm 



ionized medium (WIM) within ±100 km/s of the local standard of rest ([Haffner et al.ll2003l ) 



We used the publicly available WHAM integrated intensity data for high Galactic latitudes. The choice 
is motivated by reducing the noise from localized regions, e.g. HII regions, in the galactic disk. 

The integrated WHAM map consists of patches of contiguously observed data of size 7x7 pointings. 
The background levels in the patches are different due to different observation conditions. Such patching 
causes a systematic error, which requires a proper handling, especially in the low-signal area near the North 
Galactic Pole (NGP). 

Our data handling included the equalization of background signal level in the adjacent patches of the 
emission map. The particular algorithm is described in the Appendix [X] Figure 1 shows the original and 
the processed data. One may notice that the original data is patchy and one expects these borders of the 
patches to contribute strongly into measured spectrum. The processed data exhibits fluctuations, but do 
not exhibit rectangular domains of varying intensity. 



According to ISchlegel et al.l (|1998[ ) , dust extinction in the direction of NGP is only 3% for Ha and is 



neglected in our calculations. 

The processed data is used to estimate spectrum parameters. We do this by probing signal fluctuations 
at two different scales with subsequent fitting of the angular spectrum model. The latter has two free 
parameters, which depend on the outer scale of turbulence and amplitude of electron density spectrum. The 
results are tested for statistical consistency. 

This approach allows us to recover information about the density spectrum despite rather small dynam- 
ical range of scales characterizing the WHAM data. 

As the actual observable for Ha data is the emission measure (EM — J r? e dr ), we needed to handle the 
case of quadratic emissivity to bind the expected to be Kolmogorov spectrum of n e to the angular spectrum 
of Ha fluctuations. 



3. Theoretical expectations 

We attempt to test whether the WHAM data is consistent with the hypothesis of the Kolmogorov 
energy cascade originating at large spatial scales. To achieve this goal we, first of all, must relate the power 
spectrum of underlying density and the angular spectrum of intensity fluctuations. The latter is usually 
measured in terms of spherical harmonics. In §3.1 and §3.2 we derive the necessary expressions and make 
use of the WHAM data in §3.3. 
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3.1. Mapping of 3D power spectrum to the angular spectrum 

Let s be a homogeneous isotropic random field in a 3D space. Then its mapping to the angular 
coordinates is given by the following expression: 

/>OC 

S(9,4)= / w{r)s{r,6,4>)dr (1) 
Jo 

where w is a window function setting the layer thickness. The correspondent correlation function depends 
only on 9: 

c(6) = (s(o,o)-s*(e,o)) (2) 

After some algebra we can get: 

, , , . sin k\J r 2 — 2rr' cos 9 + r' 2 



C{9)=Att w{r)dr w{r')dr' k 2 dkF{k)—^= - (3) 

Jo Jo Jo kyr z — 2rr'cos9 + r^ 

where F is a power spectrum of s. 

Having calculated the angular power spectrum 



Ci= [ C (9) Pi(cos 8) sm9 dO 
Jo 



where 

2 



(4) 



we obtain: 

poo / poo \ 2 

Ci = 16tt 2 J k 2 dkF(k)(j ji{kr)w{r) dr J (5) 

If we choose w in a Gaussian form with the layer thickness R 

w(r) = e~«^ (6) 

we can write 

16tt 2 f°° 2 



Ci = —— u 2 duF{u/R)Qi{u) (7) 
-ft Jo 



Ql{u)= U ji(uv)e~ v dvj (8) 
After some algebra we can find the following asymptotic expression for the kernel Q: 

2h? 



«(«) - (9) 



which holds for / > 15 uniformly over u. 

The expressions [7] and [9] give us a simple mapping from 3D power spectrum F(k) to the angular spectrum 

For example, setting F to a power law with spectral index a and outer scale L 

F (fc) = ||.e-(fr) a (10) 
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we can evaluate Eq. ([7]) analytically: 



Ci = ^^ff) ^ r (11) 

22 « a -fi + HW s 



3.2. Squared Kolmogorov field 

As the Ha emissivity is quadratic with respect to the underlying physical field (electron density), we 
have to account for the related effects. 

If S is the underlying field (homogeneous and isotropic) , our emissivity can be written in the following 
form: 

s(f) = «S) + AS(r)) 2 (12) 

while (AS) =0. If J 7 is a power spectrum of AS, we can write the following expression for the emissivity 
power spectrum F, assuming that AS is a Gaussian field: 

F(k) = (2(S)) 2 T(k)+2F q {k) (13) 

where 



F q (k) = J F{k')F(k - k') dk' (14) 

Let us find an approximation for F q . 

Let us assume that J- has Kolmogorov spectral index and a low-wavenumber cutoff at kg : 

^(*) = p5r«"£ (is) 

Then, F q can be approximated using the method of matched asymptotic expansion for the asymptotic limits 
k — > and k — > oo. The asymptotic k — > oo should decrease as k~ n / 3 with a scaling factor of 2 J F(k) dk , 
which can can be illustrated as follows. As the spectral index is greater than 3, the bulk of the integrated 
value of T exists in a region near k — 0. For very large k, the auto-convolution consists of two locations of 
large integrated value, i.e. around k' = and k' — k = 0, i.e. separated by a distance k. Interpolating with 
the correct behavior in between these two asymptotic values allows F q to be approximated ai 



:rpc 
& 



F, ( k , * H ~ — '10) 

' f N 6/ll \ li / b 



\ 2Jl J 



where 



J = ^ F 2 {k)k?dk 

•M) Jo 

Ji = / T{k)k 2 dk 
■Fq Jo 



(17) 
(18) 



1 See Fig. [2] for comparison of directly calculated F q with its approximation 
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Evaluating Jo an d J\ we can write the following expression for F q : 



where 

1 t 514J r2 

F q0 = ' 13/3 , a = 0.1842 (20) 
Applying Eq. (0 to F g , after some algebra we get 



J^R 8 / 3 1 + 3.21 



-1.43, 



&) Z"/3. (i + 5 .0(^) 2A 

Finally, using Eq. (| 1 1 1) for the first term in Eq. ([7]) when i* 1 is given by Eq. (fT3")l . we have the following 
expression for the angular spectrum: 

r, - i 79 i n 2 (g) 2 ^ fl 8/3 

O/ — 1./Z-1U • /ii/3(i +0 .5(27rfl) 2 /(iL) 2 ) 4 /3 / 99 n 

„ ^2^8/3 1 l Q 11 -1.43J4; V / 

4-9 Qfi • 1 3 • . 1+3. 21e 

iU (27r /L)2/3 /n/3(i + 5.o(27ri?,) 2 /(ii) 2 ) 4/3 



3.3. Expected spectrum of fluctuations 

In this section we calculate the expectations for the angular spectrum of Ha fluctuations assuming 
that the emissivity is proportional to the squared random field having Kolmogorov spectrum (i.e. electron 
density). We shall start from Eq. (|2"2"1) . accounting for the following expression of Ha intensitjH (Smoot 
1998): 

= 0.29- f. n\dr 1 ' 

where electron density rt e is measured in cm -3 , r in pc and / in Rayleighs. 

Assuming that (n e ) is the mean electron density, R is the scale-height of the H a layer, L is the injection 
scale of the turbulence and Fq is the amplitude of the electron density spectrum one can write for the 
spherical harmonics: 

C, - 1 42 <n e > a f„fl g /» 

°' ~~ i '^ Z / 11 /3(i +0 .5(2 7 rfl) 2 /(iL) 2 ) 4 /3 

p2p8/3 1,001 -1.43J4, l Z4 J 

_i_o /II ^0 1+3. 21e 27rij. 

* (2ir/L) 2 /3 ' ; 11 /3(i +5 .o(2ir_R) 2 /(iL) 2 ) 4 /3 

where (n e ) is measured in cm -3 , i*o ~ in m -20 / 3 , R and L - in pc. 
Let us rewrite Eq. (|24|) as follows: 



r>, — 1 49 „„5/3„8/3 . <»e) flo 

O; - 1.41/ pq a iii/3(i + o. 5(27rag/0 2 ) 4 / 3 



-0.709 • v q 1 ol 



2 8/3 8/3 . 1+3. 21e 



;ii/3(l+5.0(2Trag/i) 2 ) 4 / 3 



2 Electron temperature of 8000K and 7 = 0.9 is assumed, see I Armstrong et al. I lll995h 
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where p — FqL 5 ^ 3 , q = Rq/L and a = R/Rq. For the polar region we can assume a ~ 1. 

One can see that the the angular spectrum given by Eq. (|25|) has two free parameters, p and q, which 
can be evaluated from the observatio nal data. Other param eters are taken as follows: (n e ) = 0.010cm -3 
(jCordes fc Lazioll200i . i? = 1800 pc jGaensler et al. I l2008h . 



The range of scales that is sampled with the WHAM is rather short. Therefore we shall make our 
comparison between the data and the model using a measure of the dispersion of the signal s(f/) over a 
circular area Q a (p) of radius a and centered at p: 



D(fi) = \ 



n a (p) 



s{p")dp" dp' 



(26) 



When applied to the angular spectrum given by Eq. ([25)) . it provides 

1 



D 



(27T) 



$(a, K)C(n)dK = p/i(a, a, g) +p f 2 (a,a,q) 



(27) 



where 



$(a, k) = (1 - (£(a, k)) • (j){r, k) (28) 

and <fi(a, k) — (— Ji(ok)) 2 and r is the pixel radius, which in the case of WHAM data is 0.5°. The functions 
fi in Eq. (|27|) are given by 



fi(a,a,q) = 0.256^ J 2 |U/ , (1+0i5( ^ $(o,k)k dn 

f 2 {a,a,q) = 0.70% J 2 t ii/3 ( l +5 , 0(2ffa „ /0 2 ) 4/3 ; ^(a^)AtrfAC 



(29) 



i (l+5.0(27ra( Z /i) 2 ) 

In what follows we use Eq. (|2"7| to obtain the dispersions for different values of a 



4. Probability distributions of p and q 

For the analysis below we used two sets of maps with radiH a\ = 2.27 degrees and a 2 — 8.33 degrees. 
For those values of a the dispersions calculated using the WHAM data are D\ = 7.36 • 10~ 3 ± 0.72 • 10 -3 and 
D 2 = 10.42 • 10 -3 ± 1.05 • 10 -3 , where the instrumental error of 0.04 R is accounted for. 

To estimate the consistency of our assumption of the Kolmogorov turbulence present in warm gas 
sampled by WHAM we calculate the probability distribution of the free parameters p and q. We start from 
the mutual probability of T)\ and D 2 roughly assuming that these values are independent and have normal 
PDF: 

1 / (D 1 - D x ) 2 {D 2 -D 2 f 



P D (D 1 ,D 2 )= exp - v ; 2 ^ - v ' (30) 

l-K<j D1 <j D2 \ 2a Di Za D2 J 

Solving Eq. (PZTl) in terms of p, one gets 

p = F(a,i,ati,q,Di), i = 1,2 (31) 



This numbers are taken to get maximal dynamical range from the available data 
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wherc 



F(ai,ai,q,Di) = 



(fi{a l ,a l ,q) + 4/ 2 (a» , g) ; A) 1/2 - fi(ai,ai,q) 



(32) 



Below we shall use Eq. ([31]) to derive the probability distribution P(q,p) that will be used to estimate 
the parameters p and q and their accuracy. If we write the elementary probability in the (Di, D 2 ) space as 
dP = Pd(Di, D2)dDidD2, we get the following Jacobian of the transformation: 



dS 



g^F(a 1 ,a 1 ,q, D 1 )- B ^F(a 2 ,a 2 , q, D 2 ) 



^F(ai,ai,g,Di) - J^F(a 2 ,a 2 ,q, D 2 ) 



dDidD 2 



(33) 



which results in the expression for the probability: 

-F(ai,ai,<?,L>i) - -§^F(a 2 ,a 2 ,q, D 2 ) 



Pd{D u D 2 ) 



g^F(ai, ai, q, D x )^F{a 2 ,a 2 , q, D 2 ) 
where one should substitute D t for its expression, i.e. Di = pfi(ai,oti,q) + p 2 f 2 (cii,ai,q). 



(34) 



The practical implementation of the finding probability distributions P{p) and P{q) is illustrated in 
Figures 2 and 3. The both curves were re-normalized. Initially they had encompassed the area of 0.85, 
which meant that the model spectrum was inconsistent with the observational data with the probability of 
0.15. Using the calculated PDFs, we finally have: 



q = Ro/L = Wtf, p = FqL 5 / 3 = 0.54±2$ 



(35) 



for confidence probability 0.68. Setting Rq = 1800pc (jGaensler et al. 1 120081 ). we can find the injection scale 
L and electron density spectrum amplitude Fq: 



L = 90±t° pc 

F Q = 3±1 • 10- 4 m- 2 °/ 3 
The related spectrum is shown on the Figure [5] 



(36) 
(37) 



5. Discussion 

The Big Power Law in the Sky with the spectral index coincided with the one in Kolmogorov turbulence 
was discussed in the paper of Armstrong et al. (1995). The reliable data used in their plot reflect the density 
fluctuations at scales from 10 6 to 10 13 m measured via scintillation and the electron scattering technique. 
For larger astrophysical scales the plot in Armstrong et al. (1995) contains cither upper limits or unreliable 
data, for instance, rough estimates based on turbulent velocity. 

This paper contains the first piece of evidence that the spectrum of density fluctuations at very small 
scales shown in Armstrong et al. (1995) agrees well with the spectrum of the density fluctuations measured 
at scales of 10 17 m, if we assume that the scale height of the free-free emitting layer is 1800pc. The found 
spectrum amp l itude Fq = 3 • 10~ 4 m~ 20 / 3 is within the error bar from the value 3.2 • 10~ 4 m -20 / 3 found by 



Cordes et al.. I (|1991I) from pulsar scintillations. This is a remarkable extension of the Big Power Law in the 



Sky. 
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This is suggestive of the energy being injected at scales of 40 -j- 150pc in the Galaxy and cascading up 
to very small scales. Below we discuss whether this case is a plausible one. 

Does this picture of the large scale turbulent energy cascade look reasonable? Density information alone 
cannot answer this question. Turbulence is a dynamical process, in which density can be used only as a tracer. 
More direct information is available through velocity studies. Our estimation is a lso in agreement with t he 



expected value of 100 pc associated with supernova explosions (see for instance lHaverkorn et al.l (|2008l) ). 
Studies in Chepurnov et al. (2009) of the velocity turbulence using the Velocity Coordinate Spectrum (VCS) 
provide a good fitting of the turbulence model when the injection scale is taken to be 140 ± 80 pc (VCS uses 
Fourier-transformed over vel ocity coordinate spectral data to get analytically predictable data measure, see 



Lazarian fc Pogosvanl ((2000J).). In addition, MHD simulations in Cho & Lazarian (2003) indicate that for 
subsonic turbulence the density spectrum follows well the velocity spectrum, which is Kolmogorov. This 
is however is not true for supersonic MHD turbulence, for which, according to Beresnyak, Lazarian & Cho 
(2005) the density spectrum gets shallow. 

We feel that the issue of the spectral slope does require further studies. With our limited dynamical 
range we could test the consistency of the spectral index to the Kolmogorov one. Discontinuous structures, 
e.g. ionized ridges of clouds make the spectrum more shallow. 

What is the Mach number of the free-free emitting layer in our Galaxy? This question was addressed by 
Hill et al. (2008), who compared the PDFs obtained via MHD simulations in Kowal, Lazarian & Beresnyak 
(2007) and the PDFs of WHAM data. As a result a conclusion that the sonic Mach number of turbulence 
in free-free emitting layer, which is also frequently called the Reynolds layer, is around of 2. This is close to 
the subsonic range and therefore we do not expect to see substantial deviations from the Kolmogorov scaling 
for the random density. 

All in all, the arguments above are consistent with the idea that the large scale turbulence in the 
Reynolds layer and small scale turbulence constituting the Big Power Law in Armstrong et al. (1995) are 
the parts of the universal turbulence cascade. This is also consistent with other arguments, for instance, with 
the theoretical arguments on the energy injection in turbulence. Both of the leading ideas on the injection 
of turbulent energy, i.e. via supernovae explosions and via the magnetorotational instability, inject energy 
at large scale, e.g. larger than 30 pc. This energy is bound to create a cascade in high Reynolds number 
interstellar medium. Therefore the emergence of the extended turbulent cascade is expected. 

The present study shows the consistency of the data with the Kolmogorov cascade in interstellar gas 
spanning over 10 decades. Further studies combining various data sets, including those of velocity and 
magnetic field will clarify the nature of the turbulence cascade in the Galaxy. 
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A. Removing the systematic error 



The column density WHAM map suffers from systematic error, related to uncertainty of zero levels in 
different observation blocks. It can be corrected by introducing artificial constant shifts for each observation 
block, which are set by minimization of signal differences on block borders. 

Let use denote the initial signal level difference between i-th and j-th block by Cy, if they have a 
common bordeiQ (cy = — Cji), correction shift to be found for the i-th block by d,-, and weight for adjusting 
of involvement of the border between i-th and j-th blocks by 1/wij (wij = Wji). We also define summation 
sets as f2 = {{i,j)\3cij} and Qi = {j\3cij}. 



With these definitions a function to minimize is as follows: 

I 



L = 

n y 

Taking derivative over d n , we have at the minimum: 



(°i: 



— -T—(c 



what gives us a set of linear equations for df. 



d n y o 



E C nj 
W 2 ■ 



(Al) 



(A2) 



(A3) 



To make these equations linearly independent we need to set one of dj's to some pre-defined value or 
set some other condition making the total base level of the map definite. 

In our calculations we take uiij proportional to the mean signal over the involved blocks to correct the 
algorithm's tendency to solve the problems of high-magnitude blocks at a cost of distortion of low-signal 
areas. 

Now we shall consider the calculation of Cjj. Let us denote a set of points, that belongs to the i-th block 
by A and set of points, that belongs to the j-th block by B (we take only points in a small enough vicinity 
of some point on the border, so that such set of vicinities will give us a correspondent set of estimations of 
Cij). The signal y(r) is defined on A U B, while its values on the sets A and B differ by a constant. With 
exception of this difference, we consider y{r) to be continuous and smooth enough. 

To estimate this difference on a discrete grid {r{\ we shall approximate y using a polynomial functional 
basis {/n(r)}) n — l-..iV, with account for the fact that the constant term is different in A and B (having 
the values ca and cb respectively). The other expansion coefficients c„ are equal for the both sets. 

Having applied the method of least squares, we obtain the following system of N + 2 equations for c n , 
ca and cb- 

^ ^2auB fn(l'i)fm(l'i) Y2a fm( r i) 12b fm(fi) X^AuS^/™( r *) ^ 



V Efl/n^i) 







T,BVi 



(A4) 



4 otherwise Cij is undefined 
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Having solved it, we take ca — cb as an estimation for . For the final value of Cij we take an average 
of different estimations, calculated along the border. 
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Fig. 1. — The part of WHAM map near NGP (the center of map), raw and cleaned. Map radius is 25°. 
Longitude I = points to the down side of the map. 
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Fig. 2. — F q , direct calculation (solid line) and approximation (dashed line) 
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Fig. 3. — The probability distribution function of the spectrum parameter p 
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Fig. 5. — WHAM estimation for electron density overplotted on the figure of the Big Power Law in the sky 
figure from Armstrong et al. (1995). The range of statistical errors is marked with the gray color. 



